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Abstract 

Taking into account the effect of final state interaction, we calculate the non-zero 
(naive) T-odd transverse momentum dependent distribution /ij-(x,kj_) of the pion 
in a quark-spectator-antiquark model with effective pion-quark-antiquark coupling 
as a dipole form factor. Using the model result we estimate the cos 2(fr asymmetries in 
the unpolarized ir~N Drell-Yan process which can be expressed as xhj~. We find 
that the resulting h± (x, kj) has the advantage to reproduce the asymmetry that 
agrees with the experimental data measured by NA10 Collaboration. We estimate 
the cos 20 asymmetries averaged over the kinematics of NA10 experiments for 140, 
194 and 286 GeV ir~ beam and compare them with relevant experimental data. 
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1 Introduction 



Recently it is demonstrated that the effect of final state interaction (FSI) or 
initial state interaction (ISI) can lead to significant azimuthal asymmetries in 
various high energy scattering processes involving hadrons [1,2]. Among these 
asymmetries, single spin asymmetry (SSA) in semi-inclusive deeply inelastic 
scattering (SIDIS) [1] and that in Drell-Yan processes [2] from FSI/ISI via 
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the exchange of a gluon, have been explored and are recognized as previously 
known Sivers effect [3,4]. This effect, formerly thought to be forbidden by the 
time-reversal property of QCD [5], can be survived from time-reversal invari- 
ance due to the presence of the path-ordered exponential (Wilson line) in the 
gauge-invariant definition of the transverse momentum dependent parton dis- 
tributions [6,7,8]. Along this direction some phenomenological studies [9,10,11] 
have been carried out on transverse single-spin asymmetries in SIDIS process, 
which is under investigation by current experiment [12]. Analogously the ex- 
change of a gluon can also lead to another leading twist (naive) T-odd dis- 
tribution hi(x, k^_): the covariant transversely polarization density of quarks 
inside an unpolarized hadron. This chiral-odd partner of Sivers effect func- 
tion, introduced first in Ref. [13] and is referred to as Boer- Mulders function, 
has been proposed [14] to account for the large cos 20 asymmetries in the un- 
polarized pion-nucleon Drell-Yan process that were measured more than 10 
years ago [15,16]. Recently h^(x,k\) of the proton has been computed in a 
quark-scalar diquark model [9,17] and also used to analyze the consequent 
cos 20 azimuthal asymmetries in both unpolarized ep SIDIS process [9] and 
unpolarized pp Drell-Yan process [17] respectively. 



The same mechanism producing T-odd distribution functions can be applied 
to other hadrons such as mesons. In a previous paper [18] we reported that 
non-zero h± of the quark inside the pion (denoted as h^ n ) can also arise from 
final state interaction, by applying a simple quark spectator-antiquark model. 
Among the phenomenological implications of the function hj^ n is an important 
result for the cos 20 azimuthal asymmetry in the unpolarized tt~N Drell-Yan 
process [15,16], which can be produced by the product of hy of the pion 
and that of the nucleon. Therefore, one can investigate how the theoretical 
prediction of the asymmetry is comparable with the experimental result, as 
a test of the theory and the model. In the present paper, based on from 
our model calculation, we analyze the cos 20 azimuthal asymmetry in the 
unpolarized w~N Drell-Yan process by considering the kinematical region of 
NA10 experiments [15]. To obtain the right Qt dependence of the asymmetry, 
we recalculate hj; % (x, k^_) in a spectator model similar to the model used in 
Ref. [18]. The difference is that here we treat the effective pion- quark- ant iquark 
coupling dipole form factor g n (k 2 ), in contrary to the treatment in 

Ref. [18] where we take constant. We find that h^ n (x, k^_) resulting 

from the new treatment together with h^(x,k. 2 ± ) for the nucleon in a similar 
treatment [10] can reproduce the cos 20 asymmetry which agrees with NA10 
data. We give the asymmetries predicted by our model averaged over the 
kinematics of NA10 experiments for 140, 194 and 286 GeV ir~ beam and find 
that the energy dependence of these asymmetries is not strong. 
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2 Non-zero h^. of the pion in spectator model 



In this section, we will show how to calculate h^^x, k^_) in a quark-spectator 
antiquark model. We follow Ref. [1] to work in Abelian case at first and then 
generalize the result to QCD. There are pion-quark-antiquark interaction and 
gluon-spectator antiquark interaction in the model: 

£/ = -flW^Ts^W - + h.c, (1) 



in which g w is the pion-quark-antiquark effective coupling, and e2 is the charge 
of the antiquark. When the intrinsic transverse momentum of the quark is 
taken into account, as required by T-odd distributions, the quark correlation 
function of the pion in Feynman gauge (we perform calculation in this gauge) 
is [7,8]: 



(2tt) 3 



(X) 

t+=o 



where C a (a , oo ) is the path-ordered exponential (Wilson line) accompanied 
with the quark field which has the form 



£o(0, oo) = V exp 



/ 



V 



-ig J A + (0,C,O±)df 



(3) 



etc. The Wilson line has the importance to make the definition of the distri- 
bution/correlation function gauge- invariant. Without the constraint of time- 
reversal invariance, in leading twist the quark correlation function of the pion 
can be parameterized into a set of leading twist transverse momentum depen- 
dent distribution functions as follows [13,19] 



$(x, kjj 



er l ,„kf n v 



(4) 



where n is the light-like vector with components (n + , ri~, nj_) = (1,0, 0j_), 
C/Lti/ = I[7m' 7f] anc ^ M n is the pion mass. Knowing <& n (x, kj_), one can obtain 
these distributions from equations 



f l7T (x,k 2 ± ) 
2^(x,kj)kY 

M- 



Tr[$(x,k ± ) 7 + ]; 
Tr[$(£,k ± K + ] 



(5) 
(6) 
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Fig. 1. Effective correlation function <J? in the antiquark spectator model with fi- 
nal-state interaction modelled by one gluon exchange. 

The calculation of unpolarized distribution function fi n in the antiquark spec- 
tator model can be done [20] from the lowest order (without Wilson line) cor- 
relation function in Eq. (2). However it can not lead to any T-odd distribution 
function such as hj^. As demonstrated in Ref. [1], the non-zero T-odd distri- 
bution requires final state interaction from gluon exchange between the struck 
quark and target spectator. Here we follow the observation in Refs. [6] and [21] 
that final state interaction in an initial hadron state can be taken into account 
effectively by introducing an appropriate Wilson line in the gauge-invariant 
definition of the transverse momentum dependent distribution function, or 
equivalently, quark correlation function of the hadron. The Wilson line can 
provide nontrivial phase needed for T-odd distribution functions. Since we 
have defined such a correlation function in Eq. (2), we can start from Eq. (2) 
to calculate hj; n with the explicit presence of the Wilson line. We expand the 
Wilson line to first order corresponding to one gluon exchange. Therefore, ac- 
cording to Eq. (6) and Eq. (2), can be calculated in the antiquark spectator 
model from the expression 



2/4(*,ki)ki 



x(q s \ I -ie x J A+(0,r,0±R- ] v&MOW 



h.c. 



(7) 



in which \q s ) represents the antiquark spectator state with spin s, and e\ is 
the charge of the struck quark. 

Fig. 1 is the diagram equivalent to Eq. (7). The figure shows the effective 
correlation function in the antiquark spectator model with the Wilson line 
expanding to the first order. h^ n can be obtained from the diagram by in- 
serting a l+ , according to Eq. (6). Fig. 1 is similar to the diagram used by Ji 
and Yuan [21] to calculate /fy of the proton in scalar diquark model. The 
7 5 inside the circle denotes that the pion-quark-antiquark coupling is pseu- 
doscalar coupling. The double line represents the eikonalized quark propagator 
(eikonal line), which is produced by the Wilson line along the light-cone vec- 
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tor n M = (n + ,n~,n_i_) = (0,1, 0j_). The eikonal line gives rise the final state 
interaction effect between the fast moving struck quark and the gluon field 
from target spectator system [6,21]. The Feynman rule for the eikonal line is 
i/(<? + + ie) [22] (see also Appendix in Ref. [17]), where q is the momentum 
of the gluon attached to the eikonal line. The Feynman rule for the eikonal 
line-gluon vertex is e\n^ [22]. The straight line cut by the vertical dashed line 
denotes the on-shell spectator antiquark state v s or v s . 

Usually there are two choices of the pion-quark-antiquark coupling g v : 

• case 1 : normalization constant which is used in Ref. [18]. A similar 
treatment for the proton-quark-diquark coupling g has been adopted in 
Refs. [1,17] to estimate f^(x,k 2 ± ) and /^(x, k^_) of the proton. 

• case 2 : g n as a dipole form factor [20] 

/,9n » t k 2 — m 2 , , 9 k 2 — m 2 

g n (k 2 ) = N n — — = N v (l - x" 2 



"(A 2 -P) 2 " v ; (k 2 



2\2 ' 



with 



L 2 n = (l-x)A 2 + xm 2 -x(l-x)M 2 , (9) 
k 2 ± = -(l-x)k 2 -xm 2 + x(l-x)M 2 , (10) 

and N n is the normalization constant, m is the mass of the quark/ antiquark 
inside the pion, A is the cut off parameter of the quark momentum. This 
kind of treatment has been applied to model T-even nucleon distribution 
functions [20], and recently in the calculations of frr(x, k 2 ± ) and h^(x, k^) 
of the proton [10] in order to eliminate the divergences in the /c^-moments 
of these /c^-dependent distribution functions. 

In the previous paper [18] we performed a computation on hj^ in case 1 which 
yields: 



ki) = K *, „> , (ii) 




and the corresponding unpolarized distribution is 

^(x,kl) = a ( ^+^ )2 , (12) 



where 



2 I I 

K = ^^J. mM^l-x), B w = m 2 - x{\ - x)M 2 , (13) 



n 2(2vr) 3 4vr 
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C n = (l-x)g 2 J[2(2n)% 



D n = (l+xfm 2 . 



(14) 



An interesting result is that the transverse momentum dependence of hj^.(x, k]_) 
in this model is the same as that of h± of the proton in the quark-scalar di- 
quark model [17]. 

Now we perform the computation of h\^ in case 2, that is, in the situation of g n 
as a dipole form factor. According to Eq. (7), also with the help of Fig. 1 and 
the Feynman rules introduced above, we can calculate from the integral: 



2/^(2, kj_)k 4 _L «|eie 2 | ^ f d A q _ s 2 ft + m 



M n 8(2tt)%1-x)P+^J (2 

^9n({k + q) )75- 



f d^q _ s n2s ft + m i+ 



{k + q) 2 — m 2 {k + q — P n ) 2 — m 2 + ie 

x-f + v s —^—. — ^— + h.c. (15) 
q + + ieq z — ie 

In above equation we have used (P n \ijj(0)\q s ) = v s g n {k 2 )^^i{ft + m)/{k 2 — m 2 ), 
etc, which is a result of the spectator model [20]. The 7 + in the last line 
of Eq. (15) comes from the contraction of the eikonal line-gluon vertex and 
the gluon-antiquark vertex: n il g ta/ r f = 7 + . The loop integral over the gluon 
momentum q is similar to the integral for calculating in [17] and [21]. The 
q~ integral is realized from contour method, and q + integral can be done by 
taking the imaginal part of the eikonal propagator: l/(q + + ie) —> —iii5(q + ), 
since the real part of the propagator is cancelled by the Hermitian conjugate 
term. After performing the integral we yield hj^ with a form different from 
Eq. (11): 

hiA*,*±)- 4n 2(27r)3/4(ki + L 2 )3- ^ 



The corresponding unpolarized distribution is 

fl (x k2) _ ^ 2 (l-^) 3 (kj + ^) , m 



To calculate the trace in the nominator of Eq. (15) we take the spin sum of the 
antiquark state as J2s yS v s = (tfn — ft — m), which is a little different from the 
spin sum adopted in Ref. [20]. One can find that the form of Eq. (16) is similar 
to hj; of the proton computed in Ref. [10]. We also calculate hj^., the T-odd 
distribution of the valence antiquark inside the pion, and yield h^ = hj; n . 
Comparing the two versions of h^ in Eq. (16) and Eq. (11), we find that each 
one has a significant magnitude, which means both of them can give unsup- 
pressed cos 20 asymmetry. However, the transverse momentum dependence of 
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Fig. 2. Angular definitions of the unpolarized Drell-Yan process in the lepton pair 
center of mass frame. 

the two versions are very different, that is to say, the Qt behavior of the cos 20 
asymmetry predicted by the two cases should be different. One may expect 
experiments to make a discrimination between the two versions of h^. We will 
give a further comparison with available experimental data in next section. 



3 The cos 20 asymmetry in the unpolarized ix iV Drell-Yan process 



The unpolarized Drell-Yan process cross section has been measured in muon 
pair production by pion-nucleon collision: 7r~N — > fi + fi~X, with N denoting 
a nucleon in deuterium or tungsten and a ir~ beam with energy of 140, 194, 
286 GeV [15] and 252 GeV [16]. The general form of the angular differential 
cross section for the unpolarized Drell-Yan process is 

( \ + A cos 2 9 + (j, sin 29 cos + - sin 2 9 cos 2(p] , (18) 



adtt 4tt A + 3 V 2 



where is the angle between the lepton plane and the plane of the incident 
hadrons in the lepton pair center of mass frame (see Fig. 2). The definition of 
the lepton plane depends on the choice of axes z in the lepton pair center of 
mass system. In our calculation we choose z parallel to the bisector of P^ and 
—Pat, which is referred to as Collins-Soper frame [23]. The experimental data 
show large value of v near to 30% in the Collins-Soper frame. The asymmetry 
predicted by perturbative QCD is expected to be small [15,24]. Several the- 
oretical approaches have been suggested to interpret the experimental data, 
such as high-twist effect [25,26] and factorization breaking mechanism [24]. A 
natural explanation has been proposed by Boer [14] that the product of two 
T-odd chiral-odd hy can give cos 20 asymmetry without suppression by the 
momentum of the lepton pair. In that paper, a parametrization of in a 
similar form of Collins fragmentation function [5] has been given to fit the 
experiment data. Recently it is found that non-zero hj~ can arise from final 
state interaction without violation of time-reversal invariance, and has been 
used to estimate cos 20 asymmetry in the unpolarized pp — > 11 X Drell-Yan 
process [17]. 
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Fig. 3. The cos 20 asymmetry (solid line) in the unpolarized ir~N Drell-Yan process 
defined in the Collins-Soper frame at x = x = 0.5. The dashed line represents 
asymmetry given by the previous model result [18] (with h± n calculated in case 1). 
The data are taken from Ref. [15] at 194 GeV. 

Encouraged by this proposal, we apply hf- n given by our model calculation 
to estimate the consequent cos 20 asymmetry in the unpolarized ir~N Drell- 
Yan process measured by NA10 Collaboration. In case the vector boson that 
produces the lepton pair is a virtual photon, the leading order unpolarized 
Drell-Yan cross section expressed in the Collins-Soper frame is [14] 



da{h x h 2 -> IIX) 



d^ldxidx 2 d 2 q^± 3Q 2 „ 
x cos 20 T 



J2\A(y)F[f?ffl + B(y) 



i j_,a 7 _L,a 

M X M 2 



(2h • p ± h • k ± - p ± ■ k ± 



(19) 



where Q 2 = q 2 is the invariance mass square of the lepton pair, is the 
transverse momentum of the pair, and the vector h = q±/Qt- We have used 
the notation 



FifJi] = / d 2 p ± d 2 k ± 6\ P± + k ± - qj/ifo pl)/i(£, k^ 



(20) 



From Eq. (19) one can give the expression for the asymmetry coefficient v [14]: 



(2h • p ± h • k ± - p ± ■ k ± ) 



7 _L,a 7 _L,a 

hi h-i 
M 1 M 2 



£4f[#/?]. (21] 



The cos 20 dependence as observed by the NA10 Collaboration does not show 
a strong dependence on A [15], i.e., the asymmetry is unlikely associated with 
nuclear effect. The leading contribution which comes from the valence quarks 
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Fig. 4. The cos 2(f> asymmetries in the Collins-Soper frame for pion beam with 
different energy averaged over kinematics region: x < 0.7, 4.0 < Q < 8.5 GeV 
(4.05 < Q < 8.5 GeV for the 194 GeV Data) and Q > 11 GeV. The data are taken 
from Ref. [15]. 

is hi^ u x h x therefore we can adopt the w-quark dominance, i.e., we do 
not include sea quark contribution which is expected to be small. We use 
given in Eq. (16) (case 2) to estimate the asymmetry. We also need of the 
nucleon in a similar treatment with effective coupling as a dipole form factor. 
This has been done in Ref. [10]. We use this version but only include 
contribution from the scalar diquark (/ii = /i 1 "' s , with S denoting the scalar 
diquark). There are two considerations: the first is to reduce the number of the 
parameters, and the second is because of the n-quark dominance assumption. 
Based on Eq. (21) with the denominator from the same model result, we 
give the numerical estimation of the asymmetry at x = x = 0.5 in Fig. 3 
(shown by the solid curve) with experiment data from NA10 Collaboration. 
For the parameters in the expressions of hj; n and we choose: A = 0.6 GeV, 
M v = 0.137 GeV, m = 0.1 GeV, M N = 0.94 GeV, A 5 = 0.8 GeV, and 
m N = 0.3 GeV, where M N , X s and m N are the nucleon mass, the scalar 
diquark mass and the mass of the quark inside the nucleon, respectively. For 
the coupling constant |eie2|/47r we extrapolate |eie2 |/47r — > CVa s , and take 
a s = 0.3 and CV = 4/3 which are adopted in Ref. [1]. We choose the data 
at 194 GeV of Ref. [15], since the error bars of them are smallest (the error 
in Qt is chosen to be the bin size). We find that the estimated asymmetry 
agrees with the experiment data fairly well, although our estimation is crude 
since some approximations are adopted. In contrast, as shown by the dashed 
line in Fig. 3, the Qt shape of the asymmetry produced by hj^ denoted in 
Eq. (11) (case 1 with the pion-quark-antiquark coupling as a constant) is not 
consistent with experimental data (see Ref. [18]). 

We further use given in Eq. (16) and the same parameters adopted above 
to estimate the asymmetries averaged over the kinematics of NA10 exper- 
iments. For Q T Q 2 , the momentum fractions of the quarks inside the 
pion and the nucleon satisfy the relation: xx = Q 2 /s, where i/i is the center 
of mass energy of the pion-nucleon system, for instance, for 194 GeV beam 
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y/s = ^(194 + 0.94)2 - 194 2 = 19.1 GeV, and for 140 GeV, 286 GeV beam 
y/s = 16.2 GeV, 23.2 GeV respectively [15]. The data-selecting condition of the 
NA10 experiments is: x < 0.7, 4.0 < Q < 8.5 GeV (4.05 < Q < 8.5 GeV for 
the 194 GeV Data) and Q > 11 GeV. We use above kinematical constrains to 
evaluate the averaged asymmetries for tt~ beam with different energy. In Fig. 4 
we plot the asymmetries in the Collins-Soper frame versus Qt for 140 GeV, 
194 GeV and 286 GeV beam together with the data of Ref. [15]. The estimated 
asymmetries for pion beam with different energy are still consistent with ex- 
perimental data. Our estimation shows that the energy dependence of the 
asymmetries is not strong, and this agrees with the experimental observation. 



4 Summary 

The observed large cos 20 azimuthal asymmetry in the unpolarized Drell-Yan 
process indicates a substantial non-zero value for leading twist T-odd dis- 
tribution function h±(x, k^) (which is referred to as Boer-Mulders function) 
from phenomenological aspects. Theoretically it has been demonstrated that 
of the nucleon and the meson can arise from final or initial state interac- 
tion. In this connection, we have performed a calculation of of the pion 
in a simple antiquark spectator model by taking into account final state in- 
teraction, and estimated the consequent cos 20 azimuthal asymmetry in the 
unpolarized ir~N Drell-Yan process which is then compared with experimen- 
tal data measured by NA10 Collaboration. In the calculation we adopt the 
pion-quark-antiquark effective coupling as a dipole form factor. We find that 
the resulting hj^, together with of the nucleon resulting from a similar 
treatment with nucleon- quark- diquark coupling as a dipole form factor, can 
give a good agreement of the estimated cos 20 azimuthal asymmetry with ex- 
perimental data from NA10 Collaboration. This provides a new indication on 
the role of T-odd distribution h± to the cos 20 asymmetry in the unpolarized 
Drell-Yan process from initial state interaction. 
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